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We study the occurrence of errors in a continuously decoupled two-qubit state during a \/SWAP 
quantum operation under decolierence. We consider a realization of this quantum gate based on 
the Heisenberg exchange interaction, which alone suffices for achieving universal quantum compu¬ 
tation. Furthermore, we introduce a continuous-dynamical-decoupling scheme that commutes with 
the Heisenberg Hamiltonian to protect it from the amplitude damping and dephasing errors caused 
by the system-environment interaction. We consider two error-protection settings. One protects 
the qubits from both amplitude damping and dephasing errors. The other features the amplitude 
damping as a residual error and protects the qubits from dephasing errors only. In both settings, we 
investigate the interaction of qubits with common and independent environments separately. We 
study how errors affect the entanglement and fidelity for different environmental spectral densities. 

PACS numbers: 03.67.Pp, 03.67.Lx, 03.67.-a, 03.65.Yz 


I. INTRODUCTION 

Quantum computers use superposition and entan gle - 
ment of qubits to outperform digital computers [l|, 
The advent of these machines will unquestionably en¬ 
compass a radical transformation in the way we simulate 
quantum-mechanical processes Q, imparting a plethora 
of new achievements in science and technology. However, 
to take advantage of the benefits of reliable quantum in¬ 
formation processing, we depend on the development of 
efficient ways to avoid or recover from the errors that, in¬ 
duced by environmental interaction , occur in the state 
of our quantum system. 

Accordingly, several strategies to protect quantum in¬ 
formation, particularly during the operation of a quan¬ 
tum gate, have been designed, including quantum error 
correcting codes fault-tolerant quantum computing 
[ 1 ], decoherence-free subspaces 0 , 0 , etc. One of the 
most effective methods to protect the state of a quantum 
system from decoherence is called dynamical decoupling 
(DD), which has been extensively studied in the litera¬ 
ture, both theoretically [Ml and experimentally [Ti¬ 
ll^, and is the main subject we focus in this paper. The 
DD approach is based on applying a sequence of external 
control pulses to the quantum system to be protected, 
in order to suppress the errors arising from its coupling 
with the environment. In other words, we introduce an 
additional Hamiltonian, called the control Hamiltonian, 
that acts on the Hilbert space of the system, averaging 
out the effects of the environmental perturbations. 

Alternatively, instead of control pulses, it is also 
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possible to apply continuous external fields to decou¬ 
ple the system from the environmental interactions. 
This scheme, known as continuous dynamical decoupling 
has attracted a lot of attention in recent years 
- l35| |. The CDD procedure is more experimentally 
friendly than pulsed procedures and it also sets a natural 
e for the implementation of two-qubit quantum gates 


38l| . For instance, nitrogen vacancy (NV) centers in 


diamonds have recently been shown to be strong candi¬ 
dates for use in the field of quantum technologies, with 
possible applications ranging from biological systems to 
quantum computing protocols [s^. Furthermore, for re¬ 
ducing the damage of environmental perturbations in NV 
centers, the CDD procedure has again proven to be very 
important pnl |. 


In this work, we consider CDD of a two-qubit sys¬ 
tem going through a VSWAP operation while interacting 
with a bosonic environment. While Refs, [ddl - l^ treat 
the case of magnetic noise, where the couplings between 
the qubits and their environment are orders of magni¬ 
tude stronger than the interqubit interaction, here we 
consider the case in which the two-qubit gate interac¬ 
tion is stronger than the perturbations by the environ¬ 
ment. We obtain a simple control prescription which 
allows us to prove the effectiveness of our method in a 
realistic decoherence model. The effect of the environ¬ 
ment is simulated by two different quantum channels - 
amplitude damping and dephasing - simultaneously and 
independently coupled to the qubits with different cou¬ 
pling strengths. We also study the effects of residual er¬ 
rors when the CDD protection is supplied just against the 
predominant error source, i.e., the one with the strongest 
coupling. In the present context, we consider the resid¬ 
ual error as arising from the amplitude damping chan¬ 
nel, while the qubits are dynamically decoupled from 
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dephasing. We adopt the concurrence [4T| and the fi¬ 
delity as the figures of merit of the CDD procedure. We 
show that the adopted CDD scheme provides nearly full 
protection against environmental effects when both error 
mechanisms are present and, for a residual amplitude¬ 
damping environment, a super-ohmic spectral density of 
states is more destructive than an ohmic one. Further¬ 
more, in the absence of the CDD, we see that, in the case 
of a common environment for the qubits, both entangle¬ 
ment and fidelity decay more slowly than in the case of 
independent environments. 

The present paper is organized as follows. Se ction II 
reviews the CDD procedure for the case of a VSWAP 
quantum gate. Section III describes our model to simu¬ 
late the error sources. Sect ion IV shows the results for 
the CDD protection of the VSWAP from both amplitude 
damping and dephasing, and from dephasing only, when 
amplitude damping is treated as a residual error channel. 
Finally, a conclusion is presented in Sec. V. The detailed 
solution of the master equation, that gives the dynamics 
of the reduced density matrix of the system, is presented 
in the Appendix. 


II. THE MODEL AND CONTINOUS 
DYNAMICAL DECOUPLING 

To illustrate our protective scheme, we begin by assum¬ 
ing that the interaction Hamiltonian between the qubit 
system and the rest of the universe is of the form 

= ( 1 ) 

where for s = 1,2, with Xi = x, 

(s) 

X 2 = y, X 3 = 2, and Bm , for s = 1,2 and m = 1,2,3, 
are Hermitian operators that act on the environmental 
Hilbert space. The main approach of the dynamical de¬ 
coupling method in order to reduce errors on the 

system, is to eliminate the effect of the interaction Hamil¬ 
tonian by an external control Hamiltonian. Mathemati¬ 
cally, such a condition can be written as 

Ul{t)HintUc{t)dt = 0, (2) 

Jo 

where Uc{t) is the time evolution operator associated 
with the control Hamiltonian and tc = 2 -kJu. Equa¬ 
tion (12) results from a Magnus expansion [4^ used to 
describe the total interaction-picture evolution operator, 
in the limit in which tc —>■ 0. In this limit, only the first 
term of the expansion, given by the integral in Eq. m, 
is, in general, non-zero. Thus, in this ideal circumstance, 
by imposing Eq. ([2 we are ensuring the complete elim¬ 
ination of the perturbing interactions. Here, however, 
we use Eq. ([2 as a mere guide to develop the present 
approach, for our focus is on the realistic situation of a 
finite tc- In the present paper, we will gauge numerically 
the efficacy of the resulting approximate method. 


In order to control the intensity of the exchange inter¬ 
action, possible candidates for the physical qubits should 
be, for example, properly built tunable charge qubits [d^ . 
Although, in this particular case, physical reasoning leads 
us to assume that each qubit is coupled to its own en¬ 
vironment, we shall, for the sake of completeness, also 
study the case of a common environment. For our present 
purposes, we assume that the particular form of Eq. © 
is, in the case of a common environment, given by 

H,r.t = + cr(2)) . (AH + \*B^) , (3) 

where = XB + X*B\ for s = 1, 2, A is an arbitrary 
complex three-dimensional vector, and H is a scalar oper¬ 
ator that acts on the environmental Hilbert space. How¬ 
ever, when the qubits are physically located sufficiently 
far apart, as for tunable charge qubits [i^, it is reason¬ 
able to suppose that their individual surroundings act as 
uncorrelated, independent environments. In that case, 
the particular form we assume for Eq. © is written as 

Hint = -h 

where H*^®\ for s = 1,2, acts on the environmental 
Hilbert space of the s-th qubit, and A*^®^ is an arbitrary 
complex three-dimensional vector for s = 1,2. 

In the interaction picture associated with Hc(t), the 
total Hamiltonian can be written as 

H{t) = Ho + He + U^cit)HintUc{t), (5) 

where Hq is the Hamiltonian that performs the de¬ 
sired gate operation we want to protect and He is the 
environmental Hamiltonian satisfying U\{t)HEUc(t) = 
He- We represent the environment of each qubit as 
a thermal bath of harmonic oscillators. In the case 
of a common environment for both qubits, we con¬ 
sider He = where uJk is the frequency 

of the fc-th normal mode of the common environment, 
and Cfc and ak^ are the annihilation and creation op¬ 
erators, respectively. In the case of two independent 
and identical environments, instead of the above we 

take He = I]s=i “1*^ where Uk is the fre¬ 

quency of the fc-th normal mode of the s-th qubit envi¬ 
ronment, and and are, respectively, the cor¬ 

responding annihilation and creation operators. The 
frequency Uk is the same for both independent and 
identical environments. Accordingly, we take Bm^ = 
J2k , where gk are coupling con¬ 

stants. 

In a previous work, we have shown that it is possible 
to use a continuously-applied external field to protect en¬ 
tangled states from errors caused by the unavoidable in¬ 
teractions between the qubit system and its environment 
[m . The question naturally arises as to whether it is also 
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possible to protect an entangling operation. In the fol¬ 
lowing, we show that, using the very same external-field 
configuration of Ref. [m, we can prevent errors from 
occurring during the application of a VSWAP quantum 
gate to a two -qubit s tate under decoherence. 

An ideal VSWAP gate is obtained by the Heisenberg 
coupling between two qubits, whose dynamics is governed 
by the Hamiltonian 

Hq = ( 6 ) 

where we use h = 1 throughout, J is the exchange con¬ 
stant, and, for s = 1,2, -I- 

where , and are the Pauli matrices act¬ 

ing on qubit s. Remarkably, it has been shown that 
the Heisenberg interaction alone is sufficient for universal 
quantum computation, without the need of supplemen¬ 
tary single-qubit operations Thus, a protective 

scheme for quantum gates based on this interaction, such 
as the Vs WAP operation, is of fundamental importance. 
Our results may hnd important applications in various 
experimental settings, such as double quantum dots [i^ 
and neutral atoms in optical lattices |48l |. where the ex¬ 
change interaction between two-qubits can be realized. 

It is possible to protect the considered quantum gate 
actualized by the Hamiltonian in Eq. m by the control 
Hamiltonian of the form 

Hcit) = . (7) 

In order for the evolution operator associated with the 
control Hamiltonian to satisfy Eq. (2) we must have [3l|: 

Ucit) = ( 8 ) 

since and commute, where 

= exp exp , (9) 


for s = 1, 2. 

Equations ([8]) and ([9|) imply the following external field 
configuration: 

= xnxOJ + TizUJ [z cos (uxUJt) — y sin {uxUJt)]. (10) 

Here u! = 27r/tc, nx and Uz V are non-zero integers, 
and tc is a constant. Such a field configuration is a com¬ 
bination of a static field along the x axis and a rotating 
field in the yz plane, and it is able to protect the evolu¬ 
tion described by the Hamiltonian of Eq. from the 
effects of a general class of errors. We can modify this 
field arrangement to be protective solely against a de¬ 
phasing channel by setting = 0 in Eq. m- In this 
case we only have a static field along the x axis, given 
by r2(t) = xUxUj, which is simpler than the field arrange¬ 
ment given in Eq. m- Moreover, the field is supposed 
to be spatially uniform in the neighborhood surround¬ 
ing both qubits, since it is not necessary to address each 
qubit independently. 


Because Eq. ([6l) is a scalar product, it is invariant 
under rotations and 


Ul{t)HoU,{t) = Ho. ( 11 ) 

This property of the Heisenberg interaction tremendously 
simplifies the quantum operations executed under the 
protection by the CDD, due to the fact that the in¬ 
tended gate operation remains intact under the action 
of the control fields. Without this rotational invariance, 
we would have to proceed as in Ref. [s^ and introduce 
an auxiliary rotating reference frame which complicates 
the procedure. Furthermore, this invariance has another 
important property: with the exact same field arrange¬ 
ment that preserves a quantum memory, we can also pro¬ 
tect the quantum-gate operation. Being able to protect 
the gate operation without the necessity of field recon¬ 
figuration is a tremendous simplification and certainly 
improves the prospects for experimental realization. 


III. PROTECTING THE VSWAP GATE 



FIG. 1: (Color Online) Amplitude Damping plus Dephasing: 
In figures (a) and (b) we show, respectively, the concurrence 
and fidelity for independent environments and, in figures (c) 
and (d), for the case of a common environment. The dot¬ 
ted (blue) line represents the dynamics of a VSWAP quan¬ 
tum gate with protection. The solid (red) and dashed (green) 
lines represent the dynamics without protection for ohmic and 
super-ohmic spectral densities, respectively. 


In order to illustrate our protective scheme, we take 
a Heisenberg interaction strength J = tt/8 (c.f. Eq. 
®) and a cut-off frequency for the spectral density, 
J(ui) = 7? exp(—cu/cue), given by WcT = 27r, as ex- 

plained in the Appendix, where r = 10“®s. We consider 
that the qubits interact with independent and common 
environments which are, in both cases, assumed to be at 
T = 0.2K, with a coupling constant -q = 1/20, as speci¬ 
fied in the Appendix. We investigate each environment 
model with ohmic (s=l) and super-ohmic (s=3) spectral 
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densities. For the protected cases, we set the external 
field parameters as rix = 287r/T and = IdTr/r (c.f. Eq. 
OT ). We choose the initial state of the system to be 
p(0) = |ti)(til- Note that the initial state is a product 
state and, in the ab sence of any decoherence, applica¬ 
tion of the v^SWAP gate to this initial state will create 
a maximally entangled state. 


A. Amplitude damping and dephasing errors 

In Fig. 0 we show the fidelity and concurrence of 
our two-qubit system, for the protected and unprotected 
cases, during the application of the VSWAP operation, 
when both amplitude damping and dephasing errors are 
present. We observe that, in the protected cases, the 
fidelity remains near unity during the whole time evo¬ 
lution, as opposed to the unprotected cases, where the 
fidelity decays for both ohmic and super-ohmic spectra, 
with a higher decay rate in the super-ohmic case. How¬ 
ever, since the value of the concurrence is near unity at 
the end of the time evolution of the protected cases, we 
conclude that the entangling operation is successfully car¬ 
ried out. In fact, in the protected cases shown, higher 
values of fidelity and concurrence can be obtained for 
higher values of Ux and n^. However, the same is not 
true for the unprotected gate operation. For independent 
environments, the concurrence presents a peak, but then 
decays to zero with a lower maximum value and faster 
decay rate for the super-ohmic than for the ohmic spec¬ 
trum. Finally, we see that the concurrence again presents 
a peak, followed by decay to a low but finite value, in the 
case of a common environment. 


B. Amplitude damping as a residual error 

In Fig. (I2|) we present the fidelity and concurrence 
of our two-qubit system, during the application of a 
Vs WAP gate, protected only from dephasing, while the 
amplitude damping channel is left open as a residual error 
channel. We maintain the same external field configura¬ 
tion introduced in the beginning of this section, but set 
71^ = 0, which makes our system vulnerable to the resid¬ 
ual errors. Interestingly, we see that whether our two 
qubits interact independent or collectively with the envi¬ 
ronment has very little effect on the fidelity and concur¬ 
rence. In both cases of interaction with the environment, 
the residual errors cause less damage to the system for 
the ohmic than for the super-ohmic spectral density. 

In Fig. (I3|) we show the concurrence and fidelity at 
the end of the quantum gate operation, i.e., at t = r, as 
functions of the amplitude-damping coupling constant A 
(see Eqs. and ®). For all possible cases of environ¬ 
mental spectral density and interaction, the concurrence 
and fidelity decrease linearly as functions of the coupling 
constant. While the case of a common environment with 
an ohmic spectral density presents the highest decay rate, 


the case of independent environments with ohmic spectra 
shows the lowest decay rates, as functions of A. However, 
it is important to note that, regarding these functions, 
there is little difference between the case of a common 
environment and that of two independent environments. 
A plausible argument to explain this behaviour would 
be to recall the fact that in our case qubit-qubit cou¬ 
pling is stronger than qubits’ coupling to the environ¬ 
ment. Therefore, any environmental effect on one of the 
qubits is quickly felt by the other qubit for both common 
and independent environments. 


(a) (b) 






t (x) t (x) 

FIG. 2; (Color Online) Amplitude damping as a residual er¬ 
ror channel: In figures (a) and (b) we show, respectively, the 
concurrence and fidelity for independent environments and, 
in figures (c) and (d), for the case of a common environment. 
The solid (red) and dashed (green) lines represent the dynam¬ 
ics of a Vs WAP quantum gate for ohmic and super-ohmic 
environments, respectively. 


IV. CONCLUSION 

In summary, we have presented a CDD strategy to pro¬ 
tect, from general and residual errors, a VSWAP quan¬ 
tum gate, which is an entangling operation realizable us¬ 
ing the Heisenberg exchange interation. The gate opera¬ 
tion is applied to a two-qubit system and the errors are 
introduced by amplitude-damping and dephasing chan¬ 
nels, resulting from interactions with bosonic environ¬ 
ments. We consider both common and independent en¬ 
vironments, together with ohmic and super-ohmic envi¬ 
ronmental spectral densities. We quantify the success 
of protection by looking at the fidelity and concurrence 
of the system during the time of gate operation. When 
both error mechanisms are present and no protection is 
supplied, we observe that the case of two independent 
environments with super-ohmic spectral densities is the 
most harmful. However, we have shown that our protec¬ 
tion scheme works very well against these errors, keeping 
the system at high fidelity and letting the entangling gate 
operation perform successfully. We have also considered 
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(a) (b) 



FIG. 3: (Color Online) The concurrence (a) and fidelity (b) 
at the end of the quantum gate operation, i.e., at t = r, as 
functions of the coupling constant, A, of the residual error 
channel. The solid (red) line represents the case of indepen¬ 
dent environments with ohmic spectra, the dotted (blue) line 
represents the case of a common environment with an ohmic 
spectrum, the dashed (green) line represents the case of in¬ 
dependent environments with super-ohmic spectra, and the 
purple (closely-dotted) line represents the case of a common 
environment with a super-ohmic spectrum. 


a residual error setting, where only dephasing errors are 
protected, while amplitude-damping errors are allowed 
to affect the system. In this case, we have seen that in¬ 
dependent or common environmental interactions show 
little difference and the case of a super-ohmic spectral 
density continues to be more harmful to the gate oper¬ 
ation as in the previous case. Furthermore, our CDD 
scheme uses the same external field configuration in the 
dynamic as in the static case, which is an important prop¬ 
erty for experimental applications. 

Finally, we would like to say a few words about the 
dependence of the concurrence and fidelity on the spec¬ 
tral function of separable and common environments. To 
start with, let us address the issue of separability. 

As a matter of fact, this dependence only shows up in 
the rate of change of those quantities. Separable envi¬ 
ronments seem to destroy them faster than the common 
ones. Although we are not going to provide any detailed 
explanation for this fact, we can argue that quantum co¬ 
herence properties between two qubits are more likely to 
be preserved by the presence of a common environment 
since it could mediate, at least in some cases, an indirect 
effective coupling between them. In other words, in com¬ 
mon environments there might be the possibility of exis¬ 
tence of some cooperative effect on top of the deleterious 
effects present in any coupling to general environments. 
On the other hand, separable environments would always 
act independently on each qubit giving them very little 
chance to preserve or develop any quantum mechanical 
coherence. 


As for the spectral function dependence, the stronger 
effect of superohmic environments as compared to that of 
ohmic environments can be simply understood in terms 
of time scales. It is a very simple matter for the reader to 
convince him or herself that the general spectral function 
J(uj) defined in the beginning of section [TIT] has a max¬ 
imum at ujm^ = suJc whose value is Jmlx = = 

rjujc{s/eY. Therefore, we see that for superohmic envi¬ 
ronments (with s = 3) J{uj) is peaked at Swc with its 

(3) 

maximum value given by JY ~ l.l?7Wc whereas these 
values are, respectively, ujc and ~ QArjujc for ohmic 
environments. Another important point to be observed 
is that since we are interested in protection schemes, and 
consequently testing our system for times in the interval 
0 < t < T = 271 juJc, the relevant frequency range for our 
analysis is w > 0JcI2.7t. 

Therefore we see that as time evolves from t = 0, the 
deleterious effects of the superohmic environment start 
to take place earlier than those of the ohmic environ¬ 
ment because the spectral weight of the former is more 
pronounced at a higher frequency than that of the lat¬ 
ter. In other words, it is the high frequency behavior of 
the spectral function which dominates any phenomena at 
this time scale. 

At this point one should naturally argue that our rea¬ 
soning does not make much sense for the present prob¬ 
lem because the form established for the spectral func¬ 
tion J{u}) is actually appropriate for dealing with phe¬ 
nomena at very long time scales (t ^ 1/wc) which are 
dominated by the low frequency behavior of The 

cutoff frequency ujc is a characteristic frequency of the 
environment which fixes the time scale of the problem. 
Nevertheless, we can still sustain our results if we take 
them as an indication that no matter what bath we have, 
it is its high frequency behavior that matters in this case. 
This means that if one is really interested in a more quan¬ 
titative analysis of the problem, a more detailed account 
of the environment’s high frequency behavior is in order. 
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Appendix: Solution of the master equation 

Starting from H[t), the Hamiltonian in the interaction 
picture is written as 


= E E E (a.i) 


s=l m—1 n—1 


where = cri'’^ crn\t) = 

UQ{t)an^Uo(t), for s = 1,2 and n = 1,2,3, with C/o(t) = 
exp(—iiJot)- We have used Eq. ([T]) and defined the 


operators Em\t) = U^^{t)B^UE{t)., for s = 1,2 and 
m = 1,2,3, with UE{t) = exp(—The quanti¬ 
ties Ul{t)a’^^Uc{t) = Rm,nit)a\t\ for s = 1, 2 and 

m = 1,2,3, are rotations of am , whose matrix elements, 
Rm,.n.(t), are real functions of time. We proceed as in 
Ref. and assume that the absolute temperature is 
the same in the surroundings of both qubits and these 
qubits, as well as their respective environments, are iden¬ 
tical. We then write down the master equation for the 
two-qubit reduced density matrix, in the Born ap¬ 

proximation: 


J 


dpi{t) 

dt 


2 


s,s' = l n,n' — l 


E fdt' , (A.2) 

do '' 


where we have defined the coefficients 

= E E R^At)R^',n'{t')cl:';J{t,t'), 

m—1 m' — l 

for n, n' = 1,2,3 and s, s' = 1, 2, and 

ct:2M = ErE{E^^\t)pEE^:;}{t')], 

for m, m' = 1,2,3 and s,s' = 1,2. is the 

correlation function between components m and m! of 
environmental operators calculated at the same qubit po¬ 
sition, as explained in Ref. [3l|. Here, Tr^ denotes the 
trace over the environmental degrees of freedom. The 
operators an\t), for s = 1,2 and n = 1,2,3, can be 
explicitly obtained as the components of the following 
vector relations: 

Ul{t)cr^^'>UQ{t) = + b{t)cr^‘^'> 

— c(t)((T^^) X (A.3) 

and 

Ul{t)cr’d^^Uo{t) = + b{t)(T^^'> 

— c(t)((T^^^ X (A. 4) 

where a(t) = [l-|-cos(4 Jt)]/2, b{t) = [1 —cos(4Jt)]/2, and 
c{t) = sin(4Jt)/2. The environmental density matrix. 


P£;, is taken as the one for a canonical ensemble consti¬ 
tuting a thermal bath, that is, pe = ^ exp{—PH e), where 
Z is the partition function, Z = TrE [exp(—Here, 
P = l/fcsT, ks is the Boltzmann constant, and T is the 
absolute temperature of the environment. 

We can also write the correlation function as 

= r^^’^'^TrE{Ei:Ht)pEE^^}{t')], 
where = 1 for the case of a single, common en¬ 

vironment, in which case the environmental operators 
Em\t) are independent of s, and E^®’® = Ss^s' for the 

case of two identical, uncorrelated environments. Since 
we have 

k 

and, therefore, 


k 

k 

where nu = l/[exp(/3a;fc) — 1], we obtain 


dpijt) 

dt 


2 

E 

s,s' = l ' 




2 

E 

s,s' = l ‘ 
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dt' [T^^’^'\t - t')]*m^^'^{t')]^pi{t),TZ(^'>it)] 


where 


3 3 

m—l n—1 

Ti^'^\t) = ^ \gkf rik exp(-za;fct), 

k 


and 


jds,s)^^^ = Igfcl^ (1 + Uk) expjiuJkt). 

k 

In the limit the number of environmental normal 
modes per unit frequency becomes very large, we define 
a spectral density as J{uj) = ~ ^k), with 


w G [0, oo), and interpret the summations in \t) 
and \t) as integrals over uj: 


/ dujj(uj) 


exp(—iujt) 
exp(/3w) — 1 ’ 


and 


+ r^®’®) / du}J{u!)exp{iu}t), 


for s, s' = 1,2. Here we assume a cutoff frequency Wc, 
and write J{ui) = ?7(u;®/a;®“^) exp(—w/wc) , where rj is 
a dimensionless coupling constant, and the values of s 
specify the two kinds of environmental spectral density 
we treat: ohmic (s = 1) and super-ohmic (s > 1). 
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